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The plates whose numbers are less than 2227 were taken 
before the alteration in the adjustment for tilt on 1894 
September 17. Examination of the figures in the above table 
shows that this did not sensibly alter the value of the discordance 
in question. 


A Graphical Method of Solving Kepler's Equation. 
By H. C. Plummer. 

('Communicated by U. H. Turner.) 


The prominence which Dr. T. J. J. See in a recent paper 
published in the Monthly Notices gave to the Waterston-Dubois 
method of solving Kepler's Equation makes it appear likely that 
another graphical method may have some interest. The question 
of authorship to which that paper gave rise causes hesitation in 
claiming originality; but the method which is here described has 
not, so far as I am aware, been published previously. It is pro¬ 
posed to show a way of finding an approximate solution of the 
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318 Mr. Plummer , Graphical Method lvi. 6, 

equation E—M=e sin E corresponding to a given value of M, 
e being <£ *i and *9, and also to show how the error in the 
value of E so found may be discovered very simply by graphical 
means. 



Pig. 1. 

Let a circle (centre O) of unit radius roll along a straight 
line, and let a radius OQ be divided at P in the ratio el i—e. 
Let A be the point where Q was originally on the line, M the 
foot of the ordinate through P, and E the point of contact of the 
circle at the instant. 

Then 

AM = AE—ME = arc QE-OP sin POE. 


Hence if AM=M, since the radius is of unit length, and OP=e, 
arc QE—E. But if M 7 is the point on the circle such that 
QM 7 =AM=M, M' was formerly coincident with the present 
position of M. 

It thus appears that we may take an ordinary semicircular 
protractor, marking the point P which divides the radius 
through the zero reading in the ratio e : i—e. The division will 
be facilitated if the radius has a convenient scale marked along 
it. The protractor is placed in contact with a straight ledge at 
a fixed point, so that the reading at this point is M. The instru 
ment must then be rolled along the ledge until the point P comes 
on the perpendicular to the ledge at the fixed point. The reading 
at the new point of contact is the required angle E. 

The method just described has the great advantage of 
possessing the utmost simplicity of construction, requiring no 
specially prepared curve, and being practically instantaneous in 
application. On the other hand, it is only possible to obtain very 
roughly approximate results on account of the difficulty of avoid¬ 
ing an error, amounting to a considerable fraction of a degree, in 
reading off the point of the circumference in contact with the 
ledge. It is therefore necessary to slightly modify the method 
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March 1896. of Solving Kepler's Equation. 319 

in order to get more refined approximations. For this purpose 
it must be noticed (1) that since (fig. 1) the arc QE=AE, the 
readings may he taken along, or parallel to, the base, instead of 
round the circumference ; and (2) that since OE is an ordinate, 
the abscissa of 0 may be read instead of the abscissa of E. 

Let a board be covered with millimetre paper, one set of lines 
being parallel to a straight ledge at the side of the board. Let 
the length of the base (corresponding to 180 0 ) be 450 mm. (say), 
so that 2*5 mm. correspond to a degree. Now make a semi¬ 
circle of some material, stout cardboard, for instance, whose 
curved boundary is also 450 mm., so that its radius is 

4 so 

— mm. = 143 24 mm. 

7r 

If the disc is covered with millimetre paper, with one set of lines 
perpendicular to one of the bounding radii, it is easy to find the 
point P which divides the radius in the ratio e : 1— e very 
accurately. If the semicircle is now placed with this radius 
perpendicular to the base at one extremity, and then rolled along 
the base until P has the abscissa M (the extremity of the base 
being the origin of coordinates), the abscissa of the centre is E, 
and can be read off with no greater error than o°*i, provided no 
slipping of the semicircle has been allowed to take place. 

Let E 0 be the approximation to E thus found, and let 
M 0 =E 0 — e sin E 0 . Hence AM=M-M 0 can be found at once. 

Now 

E = /'(M) =/(M 0 + AM) 

=/(M„) + AM/'(M 0 ) + +.. . 

2 

= E + _(AM ) 2 e sin E 0 I + 

0 i—e cos E 0 2 ’ (1 — ecosE n ) 2 1— ecos E 0 


Hence neglecting small quantities of order higher than the first, 


E — E 0 = aE = 


AM 

l—e cos E 0 ’ 


Now 1 —e cos E 0 can be read off at once, for referring to fig. 1, 
i—e cos E 0 = OE — OP cos POE = PM 


and it is therefore only necessary to read off the ordinate of P. 
The facility with which this correcting factor can be found seems 
to show that this method may be of considerable value in those 
cases where a second approximation by graphical means is desired. 
In practice it will probably be found best to treat the new 
approximation to E in the same way if approximations still closer 
are desired. It is, however, perhaps worthy of note that the 
calculation of the term in the above series involving a small 
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quantity of the second order is made very simple by the graphical 
method. For suppose the perpendicular at P to PE meets the 
base in T (fig. 1) ; then 

TM = PM 2 /ME = (1 - e cos E) 2 /e sin E 

and the coefficient need only be estimated roughly, since it multi¬ 
plies a very small quantity. 

The roulette traced by P is of course a prolate trochoid, to 
which PE is the normal at P, and TM is a sub-tangent. The 
instant when P falls on a given ordinate is always well marked, 
since the maximum inclination to the base of the tangent to the 
trochoid is sin- 1 *?. This fact, which can be seen very easily, 
makes it clear that the error in reading off E is sensibly the same 
for ail values of M, as well as for all values of e , for sin -I e never 
exceeds 65°. This is an important point in favour of what may 
be called the trochoid method, for under certain circumstances 
the use of the curve of sines cannot give a very good result on 
account of the approach to parallelism of the curve and the line 
x — ey=- M at the point of intersection. 

The method which has been here described is fully as general 
in its application as that of Waterston, and though I have had 
no opportunity of making an extended practical use of either 
method for the purposes of comparison, there seems to be no 
reason to anticipate any other than a successful trial for the 
trochoid method. The slipping of the disc along the base is the 
one thing to be feared, but this can be avoided by careful use 
and by giving a rough edge to the disc. Besides being equally 
well adapted to cases in which only rough approximations are 
needed, or in which the highest degree of accuracy attainable by 
graphical means is required, the trochoid method promises to be 
greatly superior when a limited use only is expected, in which 
case the trouble of constructing a reasonably accurate curve of 
sines would render the older method quite impracticable. 

Hertford College, Oxford : 

1896 March 9. 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at University of Georgia on June 5, 2015 




